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WATER  VAPOR:  ITS  STRUCTURE  AND  ELECTROMAGNETIC  ABSORPTION 


I.  INTRODUCTION 

Anomalous  infrared  (IR)  abso.  ntion  by  atmospheric  water  vapor  has  been  studied  by  the 
author1'14  and  by  other  workers15'19  to.  several  decades.  In  attributing  this  absorption  to  molecular 
clusters  in  water  vapor,  I  have  attempted  always  to  interpret  my  data  using  traditional  models  of 
cloud  physics,  nucleation  theory,  and  the  kinetic  theory  of  gases  and  liquids.  During  the  past  year, 
the  cumulative  correlation  of  data  obtained  from  IR  absorption  measurements,  mass  spectra,  electrical 
conductivity  measurements  of  moist  air,  and  thermodynamic  and  cluster  models  have  indicated  that 
the  traditional  models  are  imprecise  when  applied  to  strongly  hydrogen-bonded  substances  like  water. 

J.  G.  Wilson20  suggested  that  the  classical  Thomson  equation21,22  and  the  model  for  droplet  nucleation 
are  “.  .  .clearly  incomplete  when,  as  for  water,  strongly  polar  molecules  form  an  oriented  surface  layer. . 
Only  relatively  recently  have  detailed  treatments  of  hydrogen  bonding  effects  been  undertaken,23,6  even 
though  Thomson,  C.  T.  R.  Wilson,21,22  and  other  workers  who  formulated  the  basic  concepts  that 
have  envolved  into  today  s  cloud  microphysical,  vapor  cluster,  and  nucleation  theories  did  so  with 
absolutely  no  knowledge  of  hydrogen  bonding. 

My  investigations  to  test  the  hypothesis  that  the  IR  continum  absorption  of  water  vapor 
results  from  liquid-like,  hydrogen-bonded  clusters  of  water  molecules  (monomers)  in  the  vapor  or  in 
moist  air  have  led  to  an  apparent  confirmation  of  the  hypothesis,  and  to  quantification  of  populations 
of  neutral  water  clusters  necessary  to  account  for  phenomena  observed  by  IR  and  other  techniques. 

The  purpose  of  this  paper  is  to  present  the  main  points  of  a  new  theory  that  has  evolved 
from  my  work,  with  major  emphasis  on  physical  descriptions  of  the  populations,  sizes,  and  size 
distributions  of  the  neutral  cluster  species  that  interact  with  electromagnetic  radiation  in  the  atmosphere. 
Supporting  measurements  are  described  in  the  literature.24'27  Because  these  physical  descriptions  now 
are  available,  it  is  possible  to  extend  an  analysis  of  the  electromagnetic  absorption  and  emission  of 
vapor-phase  neutral  clusters  beyond  the  confines  of  the  IR  continum  absorption  wavelength  region 
and  into,  e.g.,  the  millimeter-wave  and  microwave  regions  where  related  spectral  effects  also  have  been 
observed.5,6  Since  these  spectral  effects  presumably  also  are  due  to  the  neutral  water  cluster  distri¬ 
butions  in  moist  or  saturated  air,  a  knowledge  of  the  cluster  species  and  their  distributions  should 
enable  spectroscopists  working  in  the  millimeter-wave  and  microwave  regions  to  predictively  model 
absorption/emission  by  these  water  species  in  those  wavelength  regions. 

The  approach  taken  in  this  paper  is  to  summarize  background  information  and  new  theory 
with  appropriate  literature  references  for  the  interested  reader.  IR  data,  however,  are  discussed  in 
considerable  detail,  and  interpretations  are  supported  by  discussions  of  the  theory.  Because  of  the 
great  sensitivity  of  cluster  populations  and  size  distributions  to  temperature  and  humidity  as  they  are 
now  understood,  earlier  estimates  of  these  parameters  from  IR  and  other  data  must,  in  effect,  be 
discarded.  Values  based  on  experimental  data  of  many  kinds  that  were  only  recently  obtained  are 
given  here.  These  values,  from  their  equations,  should  be  considered  for  use  in  subsequent  work  since 
only  in  this  way  can  users  become  familiar  with  the  extent  of  agreement  between  theory  and  their 
own  observations. 


2. 


THEORY 


2.1  General. 

The  new  theory  incorporates  several  elements  that  are  decidely  at  odds  with  much  of  the  present 
thinking  on  the  nature  of  water  vapor.  These  elements  are  first  presented  in  general  terms  and  later  defended 
by  presentations  of  the  experimental  evidence. 

Water  vapor  and  moist  air  contain  enoumous  populations  of  electrically  neutral  molecular 
clusters  which  are  present  in  peaked  statistical  size  distributions  with  typical  mean  sizes,  cu>  of  10-30 
or  more  water  molecules  (monomers)  per  cluster.12'*4  These  neutral  water  clusters  are  formed  by 
evaporation  (distillation)  of  bulk  liquid  water,  water  solutions,  or  droplets.  Only  when  liquid  water  is 
present  in  the  system  are  the  equilibria  described  in  this  paper  achieved.  When  liquid  water  is  com¬ 
pletely  removed  from  the  system,  the  neutral  cluster  populations  decay  with  half-life  times  on  the 
order  of  1 5  minutes  to  an  hour  or  more,25  and  thus  the  history  of  the  sample  becomes  important  in 
determining,  for  instance,  its  IR  spectral  properties.  Intermolecular  and  rotational  modes  of  these 
neutral  cluster  distributions  account  for  much  (probably  most)  of  the  IR  continuum  absorption  of 
water  vapor  and  can  precisely  explain  the  amplitude  and  wavelength  dependencies  of  this  absorption, 
as  will  be  shown  later,  although  direct  calculations  of  the  absorption  can  be  extremely  complex 
because  of  the  complexity  of  the  variables.  The  neutral  clusters  dissociate  to  an  extremely  small 
extent  at  any  instant  of  time13’14'25-26  to  yield  ions  by  reactions  that  might  be  represented,  e.g., 
by: 


<H2°Wl  - H+(H20)c  +  0H-(H20)d  (1) 

or 

HX(H20)c  - H+(H20)c  +  X-  (2) 

where  X*  could  be  a  negative  ion  or  an  electron.  The  energetics  of  these  reactions,  of  course,  decidedly 
favor  the  reverse  direction  as  indicated  by  my  earlier  discussion  of  expected  cluster  and  ion  mechanisms.10 
But  subsequent  experiments  with  corona  discharges,1 2  wherein  it  was  deliberately  attempted  to  enhance 
neutral  cluster  formation,  produced  no  conclusive  evidence  of  this.  Instead,  all  data  indicated  that 
the  ions  came  from  the  dissociation  of  neutral  clusters  already  present  in  the  vapor.  This  dissociation 
can  be  enhanced  by  using  an  ionization  source,  for  example,  a  /3-emitter.12  In  fact,  radioactive  dis¬ 
integration  is  the  source  of  much  ionization  in  tropospheric  moist  air.28  Thus  the  dissociation  of  the 
neutral  clusters  can  be  monitored  by  mass  spectrometry  of  ^-irradiated  water  vapor.24'26  The  electrical 
conductivity  measurements  of  water  vapor  or  moist  air  also  involve  the  mobilities  of  the  dissociative 
ioas  of  the  neutral  clusters.  Hence,  these  mobilities  also  are  a  measure  of  the  cluster  size  (c),  because 
larger  clusters  (as  ions)  move  more  slowly  in  an  electrical  field.  Measurements  recently  made29 
involving  these  mobilities  have  shown  that  it  is  very  likely  that  the  most  common  dissociative  reaction 
for  neutral  clusters,  whether  in  normal  moist  air  or  for  enhanced  dissociation  due  to  a  ^-emitter, 
simply  is: 

(H20)c+1  ■  —  —  H+(H20)c  +  OH-  (3) 

Thus  the  ion  pain  formed  by  dissociation  are  greatly  mismatched  in  size  and  mass. 


Actual  mass  spectra10,12  for  ^-irradiated  moist  air  seem  to  bear  out  this  hypothesis.  For 
example,  figure  1  shows  the  beautiful,  near-Gaussian  dissociative  ion  distributions  for  species  H+(H20)C 
from  Eq  3  for  several  partial  pressures  resulting  in  the  saturation  ratios  “s”  (=  %RH/10Q)  shown,  at 
100° C.  The  number  of  monomers  per  cluster,  “c,”  is  shown  on  the  abscissa  and  in  Eq  3.  Infrared 
measurements  show  that,  if  the  neutral  clusters  producing  the  ions  of  size  c  in  figure  1  are  taken  as 
having  c  +  1  size,  virtually  perfect  agreement  is  obtained  between  computed  spectra  and  the  IR 
continuum  absorption  spectrum.  But  the  existence  of  the  neutral  clusters  is  essential  to  account  for 
the  observed  IR  absorption  because  the  population  of  ions  is  far  too  small  to  do  so,  as  will  be  shown. 
Mass  spectral  measurements  of  the  negative  ions  are  planned.  Note  the  higher-than-expected  “spike” 
at  c  =  21  in  the  upper  spectrum  of  figure  1.  This  corresponds  to  the  population  of  neutral  cluster 
having  the  favored  low-energy  configuration  of  a  central  monomer  surrounded  by  a  symmetrical 
pentagonal  dodecahedron  23  which  begins  to  approach  a  spherical  shape.  Such  spikes  always  are 
observed  in  these  spectra  immediately  to  the  left  of  a  lower-than-expected  population,  e.g.,  c  =  22  in 
the  upper  spectrum  of  figure  1.  This  allows  the  interpretation  that  the  neutral  clusters,  except  for 
very  small  open-chained  ones,  always  have  configurations  such  that  one  or  more  favored  geometric 
solid-like  structures  are  chained  together  or  have  tails  of  OH'  groups,  or  longer  chains,  attached  to  a 
hydrogen  atom  imbedded  in  the  cluster  cage  matrix.  Thus  c  =  21  is  more  populous  because  it 
corresponds  to  an  OH'-depleted  population  added  to  the  expected  c  =  21  cagelike  population  already 
present.  This  is  exactly  the  sense  of  Eq  3. 

Further  defense  of  this  interpretation  of  mass  spectra  is  as  follows.  The  technique  of 
studying  vapor-phase  equilibria  by  the  dynamic  method  of  rapidly  cooling  gases  to  “freeze”  the 
equilibrium,30  that  is,  to  prevent  the  shift  of  equilibrium,  has  been  accepted  in  physical  chemistry 
since  at  least  as  early  as  1910  when  Koref31  reported  its  use  to  obtain  very  precise  equilibrium 
constants  for  carbon  disulfide.  Precisely  this  kind  of  “freezing”  takes  place  when  a  water  vapor  sample 
is  subjected  to  rapid  expansion  and  cooling  from  near-atmospheric  pressure  through  a  tiny  orifice  and 
into  the  high  vacuum  typical  of  mass  spectrometers.10  But  it  has  been  argued  that  mass  spectra  of 
water  clusters,  although  they  might  contain  related  information,  cannot  be  taken  as  direct  representations 
of  cluster  ion  distributions  in  moist  air  samples  from  which  they  are  obtained,32'33  primarily  because 
the  stresses  of  rapid  expansion  and  extreme  cooling  (freezing)  distort  the  resulting  distributions  that 
are  observed.  Certainly,  it  would  be  difficult  to  rationalize  mass  spectra  obtained  for  typical  moist 
air  samples  if  it  were  assumed  that  they  must  form  by  the  mechanism  of  monomer  collision  according 
to  classical  kinetic  theory  of  gases  and  were  then  ionized  by  a  radiation  source  in  the  mass  spectrome¬ 
ter.  The  hydration  of  ions  to  form  heterogeneous  clusters  that  can  be  accurately  studied  by  mass 
spectrometry  is  well  documented,34  but  the  situation  is  much  less  clear  for  plain  water  vapor  where 
it  must  be  assumed  by  kineticists  that  water  monomers  first  are  ionized  and  then  hydrated  to  produce 
the  observed  mass  (size)  spectra.  Indeed,  a  much  simpler  evaporative  mechanism  to  explain  the  cluster 
distributions  in  the  vapor,  and  their  resulting  mass  spectra,  can  be  shown  to  produce  excellent 
agreement  between  theory  and  experiment. 

Thus  the  dissociative  ion  cluster  mass  spectra  of  figure  1  can  be  argued  to  be  nearly 
identical  to  the  size  spectra  of  the  neutral  cluster  populations  that  produced  them.  If  Eq  3  were 
correct,  the  spectra  of  figure  1  would  be  precisely  those  of  the  neutral  clusters  if  the  abscissa  were 
relabeled  (c  +  1)  instead  of  (c).  This  would  explain  the  excellent  agreement  between  the  dissociative 
ion  mass  spectra  and,  for  example,  IR  absorption  trends  with  temperature  and  humidity.10'12 
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Figure  1.  Mass  Spectra  of  ^Irradiated  Moist  Air  at  Partial  Water  Vapor  Pressures, 
p,  mm  Hg  (Torr)  and  Saturation  Ratios,  s  =  %RH/ 100,  Shown 
All  spectra  are  for  100°C,  showing  relative  populations  of  H+(H20)C  species  (Eq  3). 


From  Eq  3,  the  dissociative  equilibrium  product  of  the  neutral  clusters  can  be  written: 

a  H+(H20)c,aOH- 

K  =  -  (4> 

a(H20)c+1 


where  a  is  the  molar  concentration  of  any  species.  For  water,  either  in  the  liquid  or  the  vapor  phase, 
the  magnitude  of  the  denominator  in  Eq  4  remains  almost  constant  because  of  the  small  extent  to 
which  water  dissociates  compared  to  other  liquids  or  solutions.  Thus  the  so-called  ion  product,  Kw, 
of  water  is  used  instead  as  the  equilibrium  dissociation  constant,  and  it  is  widely  tabulated.  This 
product  is  also  capable  of  being  used  in  equations  to  model  the  temperature  dependence  of  the  IR 
continuum  absorption. 8,9,1 8  The  ion  product  is: 


Kw  =  aH+(H20)c‘a°H'  (5) 

If  we  now  designate  the  total  population  or  number  of  neutral  water  clusters  of  all  sizes  in  1  cc  of  water 
vapor  or  moist  air  as  Ncc,  and  the  number  of  dissociative  ions  of  these  clusters  per  cc  at  any  instant  as 
Icc,  it  is  easy  to  see  that: 

Icc  =  (k)2v/Kw(Ncc)  (6) 


where  k  is  a  normalizing  constant  and  “2” in  Eq  6  indicates  that  ion  pairs  are  formed  according  to  Eq  3. 

But  when  actual  measurement  data  are  used  to  evaluate  k  in  Eq  6,  it  is  found  that  k  ~0.5  or,  in  other  words, 
that  the  dissociative  ion  cluster  mobilities  are  greatly  different,  thus  supporting  the  earlier  discussion.  In 
fact,  it  is  the  OH'  ions  (Eq  3)  that  are  small,  fast,  and  mobile,  offsetting  the  mobilities  of  the  comparatively 
large,  slow,  and  immobile  H+(H20)c  ions  in  determining  the  measured  ion  content  of  a  vapor  sample  by 
electrical  conductivity.  Actual  measured  values  of  Icc  are  shown  in  figure  2  vs  Celsius  temperature  on  the 
abscissa.29  The  upper  curve  shows  the  total  ion  population  measured  for  fully  saturated  vapor  (  s  =  1), 
while  the  lower  diagonal  curves  show  measured  values  of  Icc  for  vapor  samples  of  constant  partial  pressure 
at  several  temperatures. 

Data  like  those  in  figures  1  and  2  can  be  combined  with  Eq  6  to  yield  directly  the  total 
population  of  neutral  clusters  per  cubic  centimeter  of  vapor,  Ncc,  vs  saturation  ratio  (=  %RH/100)  and 
temperature  for  vapor  in  contact  with  liquid  water  at  equilibrium,  as  shown  in  figure  3.  This  figure 
shows  that  at  any  constant  temperature  (diagonals),  Ncc  has  an  (s)2  dependency,  i.e.,  a  pressure-squared 
dependency  just  like  that  of  the  IR  continuum  absorption.13  Figure  3  also  shows  that  when  a  vapor 
sample  is  heated  at  constant  partial  pressure  and  progressively  lower  saturation  ratios  result,  the  value 
of  Ncc  remains  virtually  unchanged,  as  shown  by  the  hozizontal  line  at  Ncc  =  10* 5.  This  indicates 
that  Ncc  is  a  constant  for  a  given  partial  pressure  of  water  vapor  at  1  atm  total  pressure,  at  least  for 
conditions  to  the  left  of  the  nearly  vertical  diagonal  in  figure  3,  which  terminates  at  the  boiling  point. 
The  significance  of  this  nearly  vertical  diagonal  will  be  discussed  in  the  appendix. 

Clearly,  the  IR,  or  other  electromagnetic  absorption  due  to  molecular  clusters  of  water  in 
the  atmosphere,  will  be  proportional  to  the  total  cluster  population  per  cubic  centimeter,  Ncc,  times 
the  mean  size  of  the  neutral  cluster  distribution,  cu.10  Ncc  is  given  directly  for  a  wide  range  of 
conditions  in  figure  3.  Values  of  c  and  cu  have  been  measured  for  many  temperatures  and  partial 
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Figure  2.  Measured  Values  of  the  Total  Population  per  cc,  1^,  of  Dissociative  Ions 
of  Neutral  Water  Clusters  After  Eq  3 

The  circular  and  triangular  points  are  taken  directly  from  electrical  conductivity  data, 
while  the  diamond  and  square  points  are  obtained  by  calculation  from  mass  spectral 
data.  Ion  mobility  was  taken  as  u  =  0.6  cm^/V-sec,  and  data  are  for  I  atm  total  pressure 
in  moist  air.  Three  constant-partial-pressure  curves  are  shown,  having  shallower  slopes 
than  the  saturation  (s  =  1. 0)  curve.  These  three  curves  correspond  to  heating  saturated 
moist  air  (s  =  1.0)  to  the  temperatures  shown  from  a  starting  temperature  of  (top  to 
bottom)  54° C,  35°C  and  28°C. 


Figure  3.  Total  Population  per  cc  of  Neutral  Water  Clusters,  N^.,  Obtained  from 
Experimental  Data  for  Ion  Population,  Equilibrium  Constants  and  IR  Absorption 

Computed  using  Eq  6  from  data  like  those  in  figures  1  and  2;  the  horizontal  line 
represents  constant  partial  pressure,  while  the  nearly  vertical  curve  to  the  right  is 
explained  in  the  appendix. 


pressures  by  taking  mass  spectra10  of  moist  air  like  those  shown  in  figure  1.  Thus  if  the  total  fraction 
of  all  water  vapor  that  is  clustered  is  designated  as  (nc)v,  it  is  straightforward  that: 


(nc)y 


(NccKch,) 

N1 


(7) 


where  Nj  is  the  population  of  monomers  per  cubic  centimeter  and  is  given  by: 

<s-P°>  NAvo 

N. -  (8) 

1  R0 

where  R  is  the  gas  constant  =  62,360  cc-Torr/°K  g-mole  (a  Torr  is  a  mm  Hg),  0  is  the  Kelvin 
temperature,  (s  p°)  is  the  saturation  ratio  times  the  saturation  vapor  pressure  (Torr)  at  this  temperature, 
and,  hence,  the  partial  water  vapor  pressure  (torr)  and  Navo  **  Avogadro’s  number  =  6.02  X  1023 
molecules  per  g-mole. 

2.2  Data  Interpretation  from  the  New  Theory. 

Nonlinear  data  from  many  kinds  of  measurements  can  be  explained  by  the  new  theory. 

Some  examples  include: 


2. 2. 1  Infrared  absorption  measurements  produce  nonlinear  results  at  higher  saturations  such  as  those 
seen  by  Bignell. 1 5  Spectral  absorption  data  both  in  the  IR  and  at  millimeter-  and  micro-wavelengths5,6 
are  extremely  hard  to  reproduce,  especially  in  humid  climates.  Measurements  in  steam16,17  show 
different  absorption  behavior  of  water  vapor  than  that  seen  at  lower  humidities,  with  sudden  changes 
in  spectral  shapes,  luminescence-like  activity,  and  evidence  that  a  very  large  fraction  of  the  vapor, 
perhaps  several  percent  or  more,  is  clustered.7  Vapor  spectra  taken  with  droplets  present  (thus 
assuring  intimate  vapor/liquid  interaction  with  large  interfacial  areas  to  allow  continuous  evaporation 
and  recondensation)  show  intense  spectral  activity  that  could  be  attributed  to  neutral  clusters, 
especially  in  the  millimeter-wave  region18  where  droplet  scattering  is  not  confused  with  absorption. 

2.2.2  In-house  measurements  of  mass  spectra  of  moist  air10  showed  that  samples  for  analysis  at 
99-100°C  could  only  be  obtained  by  drawing  saturated  water  vapor  into  the  mass  spectrometer  from 
flasks  heated  to  a  maximum  temperature  of  85-90°C.  Attempts  to  use  warmer  saturated  vapor,  or  to 
directly  approach  saturation  at  99-100°C  in  the  mass  spectrometer  sample  volume,  produced  spectra 
for  which  the  larger  cluster  sizes  in  the  distribution  (c  =  45-50)  no  longer  gave  smooth,  Gaussian-like 
distributions  but  very  ragged  ones.  This  allowed  the  interpretation  that  nucleation  of  the  clusters  of 
sizes  45-50  was  occurring  near  the  boiling  point,  a  situation  not  allowed  for  by  classical  nucleation 
theory20  which  requires  a  so-called  “critical  supersaturation,”  (s)^-  of  — 1.9  for  neutral  cluster 
nucleation  to  form  droplets  near  100°C. 


2.2.3  My  measurements  of  Icc  (figure  2)  showed  that  enormous  total  ion  populations  were 
possible,  at  temperatures  approaching  100°C,  of  nearly  1012  per  cubic  centimeter.29  These  ion 
populations  were  sufficient  (Eq  6)  to  account  for  all  monomers  at  the  boiling  point  and  1  atm  total 
atmospheric  pressure.  Furthermore,  when  ion  populations  were  measured  at  ambient  temperature  and 
lower  saturation  ratios,  they  increased  linearly  with  s,  i.e.,  with  an  s1  power  dependency,  until  a 
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saturation  ratio  of  0.S-0.6  (50-60%  RH)  was  readied.  At  this  point,  Icc  values  measured  vs  s  began  to 
increase  markedly,  producing  a  sharp  knee  in  the  data  curves  that  gave  the  appearance  that  1^.  had 
encountered  some  sort  of  diagonal  barrier  and  was  attempting  to  climb  along  it  as  s  was  increased 
further.29 

2.2.4  The  kinetic  theory  of  evaporation  and  recondensation35  for  water  at  the  vapor /liquid 
interface  requires,  to  account  for  equilibrium  evaporation  rates  actually  measured  for  water,  a  saturation 
ratio  of  s  =  0.99999+  at  the  evaporative  boundary  layer  under  virtually  all  conditions.  This  implies 
that  at  the  boundary  layer,  the  vapor  and  liquid  phases  are  not  clearly  distinguishable  and  that,  since 
liquid  water  is  nearly  completely  hydrogen-bonded?6  the  vapor  could  be  also.  This  could  also  account 
for  the  experimental  observation  that  the  dissociative  ion  product  of  water  (Eq  5)  applies  to  either 
the  vapor  or  liquid  phase.  Croxton37  states  that  no  statistical-mechanical  treatment  of  the  vapor/liquid 
interface  is  yet  reliable,  and  that  water  presents  a  particularly  formidable  problem  in  its  interfacial 
properties. 

The  mean  size,  cu,  of  the  neutral  cluster  distribution  increases  with  increasing  s  at  constant 
temperature,10  but  under  conditions  of  constant  partial  pressure,  cu  can  be  described  as  a  function 
of  temperature  alone.  From  experimental  mass  spectra,  the  equation  of  c„  is  (see  appendix): 


(cu)s,0  -  (s-pfl)1/2exp  (1919/0  -  4.655), 


(9) 


while  the  general  equation  for  the  ratio  of  cu’s  for  two  conditions  is: 


^^2 

(®u>l 


•  1/2 
,m\  e 
ViPi  / 


exp  r  (1/^2  •  l/0j) 


(10) 


where  in  this  case  r  =  1919. 


Consequently,  the  negative  temperature  dependency  of  the  IR  continuum  absorption,  e.g.,  could  be 
determined  by  the  negative  temperature  dependency  of  the  mean  cluster  size  at  constant  partial 
pressure.  If  one  wishes  to  assume  a  Boltzmann-like  temperature  dependency  for  cu,  hydrogen  bond 
energies  on  the  order  of  -0.1  to  -0.2  ev  from  Eq  9  will  be  computed.  Results  of  such  magnitude 
often  have  been  reported  by  spectroscopists  who  assumed  that  the  dimer  (cluster  of  size  c  =  2)  alone 
accounted  for  anomalous  absorption  by  atmospheric  water  vapor,6  but  they  assumed  that  it  was  the 
population  of  dimers,  not  the  cluster  size,  that  had  this  dependency. 

From  Eq  7,  the  total  cluster  fraction  (n^y  equals  the  product  of  (N^)  times  cu  divided  by 
the  monomer  population,  N|.  Actual  values  and  ratios  of  Cy  are  given  in  Eqs  9  and  10.  From 
experimental  measurements,  the  equations  of  Ncc  also  are  known  (figure  3): 

(Ncc)s,0  »  1.056  X  10,2(s-pJ)2  01) 
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and  thus  the  general  equation  for  the  ratio  of  N^’s  for  two  conditions  is: 

*Ncc>2  /S2P2\2 
(NcPl  \81.P°1  > 


which  is  constant  at  constant  partial  pressure.  By  combining  Eqs  10  and  12,  the  general  equation  for 
the  ratio  of  the  product  (NccKcy)  for  two  conditions  is  obtained: 


<Ncc>2^)2 

(Ncch(cuh 


exp  r  (1/02  '  1/0 1> 


(13) 


From  Eq  7,  the  general  equation  for  the  ratio  of  total  cluster  fraction  for  two  conditions  is: 


(n<Pv2 

(nc)vi 


02 

—  exp  t  (1/02  *  l/0i> 
01 


(14) 


3.  THE  INFRARED  CONTINUUM  ABSORPTION 

The  neutral  cluster  species  for  which  detailed  physical  descriptions  have  been  given  earlier 
would  always  exist  in  the  atmosphere,  and  when  in  equilibrium  with  liquid  water  they  would  be 
found  in  the  populations  and  sizes  given  in  figures  1  and  3.  These  neutral  clusters  would  interact 
with  electromagnetic  radiation  at  all  wavelengths,  to  a  greater  or  lesser  extent,  owing  to  the  different 
kinds  of  cluster  modes  in  different  spectral  regions.  While  the  following  discussion  will  focus  on  the 
IR  continuum  absorption  because  the  IR  is  the  spectral  region  best  known  to  me,  a  treatment 
similar  to  the  following  could  be  used  by  spectroscopists  in  other  wavelength  regions  to  model  cluster 
activity  in  the  atmosphere. 

Since  the  electromagnetic  absorption  at  any  wavelength  due  to  neutral  clusters  in  the  vapor 
will  be  directly  proportional  to  the  product  of  cluster  population,  Ncc,  times  average  size,  Cy,  (Eq  7) 
which  is  proportional  to  the  fraction  of  the  vapor  clustered,  (nc)v,  it  follows  that: 


“v 2  ^nPv2 

®V|  ^nPvj 


(15) 


where  is  a  mass  absorption  coefficient  having  units  like  cm^/g  or  m^g.  Thus,  the  ratio  of  or’s  in  Eq  15 
also  is  equal  to  the  right-hand  side  of  Eq  14.  Furthermore,  the  so-called  self-broadening  coefficient,  (Cp, 
of  the  IR  continuum  absorption  is  related  to  a  by: 


(Cs)2  “v2/«lPl 
<cs>j  “vj  V*2P2 


(16) 


16 


where  (Cs)  has  the  units  cm2/monomer-atm.  Thus  combining  Eqs  14,  IS,  and  16  gives,  at  some 
specified  wavelength  X  (pm): 


Eq  17  is  similar  to  one  proposed  by  Roberts  etal., 19,1 3  to  model  the  temperature  dependency 

O  O 

observed  experimentally  for  (Cs)^  when  the  reference  temperature  0  j  =  296  K: 


<<>X,8  =  «$X,296  “P  ro  <'»2  -  >«««) 


(18) 


However,  Eq  18  omits  all  of  the  pressure  and  temperature  terms  to  the  left  of  the  exp  term  in  Eq  17. 
This  suggests  that  Eq  18  is  correct  only  under  certain  circumstances,  i.e.,  when  near-constant  partial 
pressure  is  maintained  between  two  temperatures  with  only  a  slight  pressure  variation  needed  to  offset 
the  ratio  of  absolute  temperatures.  We  can  examine  this  by  plotting  Eqs  17  and  18  together,  for 
r  -  tq  =  1919,  as  is  done  in  figure  4.  Eq  17  produces  the  field  of  curves  shown  when  referenced  to 
23°C  (Roberts  et  al.,  296°K).  The  reference  point  in  figure  4  is  shown  at  s  -  0.457,  which  is  on  the 
right-hand  diagonal  discussed  in  the  appendix.  Figure  4  includes  dashed  curves  showing  Eq  18  and  also 
constant  partial  pressure  conditions.  The  interpretation  is  borne  out  that  Eq  18  works  only  for  ~ 
constant  partial  pressure  (probably  due  to  the  manner  in  which  experiments  were  run  from  which  it 
was  derived),  while  Eq  17  is  general  and  should  work  very  well  for  all  combinations  of  water  vapor 
partial  pressure  and  temperature  in  the  real  atmosphere. 


In  an  earlier  paper,13  I  gave  an  equation  that  can  model  the  behavior  of  figure  4  because  it 
is  based  on  the  dissociative  ion  product  of  water  (Eq  S): 


=  (Cs}X,0j 


exp[3578  (l/0j  -  1/02)1 


(19) 


Now  it  is  possible  to  solve  for  the  exponents  m  and  n  in  Eq  19  simply  by  setting  this  equation  equal 
to  Eq  17,  from  which,  exactly: 


m  =  -  1/2 


and 


n  =  1/2 


and  thus  the  general  equation  for  figure  4  also  can  be  given  by: 


yw —  ^ 


which  for  constant  partial  pressure  gives: 


Constant  Partial  Pressure 


P°1 

T*  exp  3578(1/5,  -  l/07) 
P  2 


(21) 


Eq  21  justifies  exactly  the  original  observation8,9,13  that  (Cp^  is  directly  proportional  to  the  ion 
product  as  (Eqs  5  and  6)  divided  by  pressure: 


«?x- 


(22) 


where  pw  =  p0^  and  the  ratio  of  (Cp’s  is  given  in  Eq  21  since 


13 


, —  /-3578  -  4.085\ 

vKw  =  exP\”T  / 


(23) 


For  comparison  with  Eq  18,  which  is  referenced  to  296°K,  Eq  20  can  be  written: 


1/2 

exp  (14  -  3578/8) 


(24) 


in  general  form;  or  for  constant  partial  pressure  in  the  form: 

Constant  Partial  Pressure,  (Cpg  =  exp  (15. 14  -  3578/0)  296°K:  (25) 

Eqs  24  and  25  will  reproduce  the  curves  of  figure  4  quite  precisely.  All  of  the  equations  given  in  this 
section  can,  of  course,  be  written  to  reference  any  temperature  from  -16°C  to  100°C,  the  limits  of 
vapor/liquid  equilibria  at  1  atm. 

The  diagonal  curve  shown  in  figure  4  and  discussed  in  the  appendix  has  been  mentioned 
previously.  If  one  humidifies  a  sample  and  thus  increases  s,  for  example  at  constant  temperature,  it 
can  be  seen  that  when  the  diagonal  is  reached  the  measured  values  of  (Cp  vs  s  will  turn  sharply 
upward.  This  will  lead  the  spectroscopist  to  believe  that  he  has  caused  condensation  to  occur  on  his 
optics  or,  if  the  optics  are  heated,  that  droplet  growth  has  begun  on  conventional  condensation 
nuclei,  making  his  data  useless.  In  reality  what  has  happened  is  that  the  nuetral  cluster  population 
(figure  3)  increases  more  rapidly  with  s  under  these  conditions.  The  data  obtained  by  the  spectroscopist 
will  be  perfectly  legitimate  cluster  absorption  data,  but  will  be  confused  with  conventional  interpre¬ 
tations  that  droplets  are  present  and  are  scattering  or  absorbing  radiation.  But  the  clusters  themselves 
will  grow  into  droplets  if  allowed  to  nucleate.  The  existence  of  clusters  and  droplets  together  is  not 
only  explicable  by  the  new  theory,  it  is  indeed  inescapable. 
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Regardless  of  wavelength,  the  anomalous  absorption  of  water  vapor  can  be  attributed  to  modes 
of  hydrogen-bonded  molecular  clusters  like  those  found  in  the  liquid  phase  which  are  present  as 
fractions,  (nc)v,  of  the  vapor  given  in  several  equations  in  this  paper.  But  only  some  fraction,  which 
can  be  called  (ng)x,  of  (nc)v,  will  be  active  in  the  absorption  at  some  wavelength  X.  Spectra  of 
steam17  give  clues  about  the  wavelength-dependency  of  absorption  produced  by  the  neutral  cluster 
distributions  in  vapor  like  those  in  figure  1.  For  example,  at  400°K  data  like  those  in  figure  5  are 
obtained.  Earlier  analyses  have  shown  that  for  a  variety  of  cluster  configurations  some  modes  of  the 
clusters  should  be  resonant  at  wavelengths  such  that:8,9 


X  =  6.4  Vc  (26) 

This  seems  to  be  borne  out  in  figure  5,  where  the  four  vertical  lines  correspond  to  wavelengths, 
reading  from  left  to  right,  of  cluster  sizes  in  Eq  26  of  c  =  3,  4,  5,  6.  But  note  also  that  the  amplitudes  of  the 
recognizable  peaks,  which  are  spaced  properly  only  for  a  constant  of  6.4  in  Eq  26,  are  only  a  small  fraction 
of  the  amplitude  of  the  overall  envelope  of  the  IR  continuum  absorption.  This  suggests  that  the  modes 
responsible  for  the  relationship  in  Eq  26  are  only  some  of  those  responsible  for  the  overall  absorption  or,  at 
the  very  least,  that  the  individual  cluster  resonances  give  rise  to  absorption  bands  that  are  so  broad  that  they 
overlap  those  of  many  neighbors  and,  therefore,  that  the  continuum  absorption  envelope  adds  together  the 
“wings”  of  many  bands  or  lines.  Eq  26  turns  out  in  practice  to  be  very  useful,  and  it  can  be  easily  shown 
that  spectra  computed  by  it  from  cluster  distributions  match  the  continuum.  The  overall  absorption  will 
be  proportional  to  the  product  of  Ncc,  cu,  and  an  average  absorption  coefficient  due  to  hydrogen-bonded 
monomers  in  the  vapor  that  can  be  designated  («h^X  an<1  can  8*ven  the  units  cm^/clustered  molecule. 
There  is  no  reason  to  believe  that  this  absorption  coefficient  cannot  be  calculated  from  the  known  spectral 
properties  of  liquid  water  at  wavelength  X,  because  the  clustered  molecules  are  like  those  found  in 
liquid  water.  If  one  treats  the  clusters  in  vapor  as  tiny  aerosols,  it  follows  that  their  average 
absorption  coefficient  can  be  computed  from:11 

4tt  X  104kxflmx)M 

-  —  '•  -  -  (27) 

Xp  (nc)L  NAvq 


where  kx  is  the  imaginary  part  of  the  complex  index  of  refraction  (n  -  ik)x  =  mx,  for  water  at 
wavelength  X,  M  =  18  is  the  molecular  weight,  Navo  1s  Avogadro’s  number,  p  ~1.0  is  the  liquid 
water  mass  density  (g/cc),  (nc)L  is  the  fraction  of  liquid  water  clustered  at  the  temperature  of  interest 
( — 0. 85  at  20°C)36j  an(j  fi(mX)  is  given  by: 


f(mx)  = 


9n 

(n2  +  k2)2  +  4(n2  .  X2)  +  4 


(28) 


where  nx  is  the  real  part  of  the  complex  index  of  refraction.  Very  precise  values  of  nx  and  kx  have 
been  given  by  Hale  and  Querry.38  For  many  approximation  calculations,  one  can  take  f(mx)  ~1.0. 
From  the  above  equations  and  assumptions  it  follows  that  the  absorption  coefficient  of  water  vapor 
that  can  be  attributed  to  neutral  clusters  (Bignell’s  ‘l^”15),  (cm^/g)  iS: 

4jt  X  104(ng)x(s-p°)kxfl(mx) 

(av)X,fl  =  X-p-  (1)  (760)  (nc)L  (29 
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Hence,  the  “self-broadening”  coefficient  of  the  1R  continuum  absorption  is: 


M  (1X760) 

lC?x  '  * 


104(ng)x 


M  k ^  ffm^) 
NAvo  x*>  (nc>L 


(30) 


Roberts  ef  al.  19  have  given  a  useful  equation  for  the  spectral  dependency  of  the  IR 
continuum  absorption  which  can  be  rewritten  for  wavelength  instead  of  wavenumber  as: 


(C?X,296 


*  1 


=  10‘22  1.25  +  1670 


(31) 


for  0  =  296°K.  Using  Eq  31  with  Eq  30  and  taking  p  =  f(mX)  =  1.0,  one  can  calculate  the  fraction 
of  (nc)v  that  is  active  in  the  absorption  at  wavelength  X,  i.e.,  (ng)^.  The  result  for  296°K  is  shown  in 
figure  6.  The  results  are  consistent  with  the  earlier  observation  that,  even  in  stream,  the  fraction  of 
IR  absorption/emission  attributable  to  cluster  modes  active  near  X  =  10  um  is  only  a  few  percent,7 
even  though  nearly  all  water  vapor  is  clustered  near  the  boiling  point.  Near  the  peak  of  the  IR 
continuum  absorption  beyond  X  =  30  pm,  figure  6  indicates  that  almost  all  clustered  molecules 
contribute  to  this  absorption. 


But  there  is  no  need  to  show  that  the  entire  continuum  absorption  through  the  rotational 
band  at  30-50  pm  is  due  to  clusters.  Much  of  this  absorption,  especially  at  longer  wavelengths,  can 
be  attributed  to  monomer  modes.  Thus  there  are  alternative  explanations  ot  the  ranging  of  (ng)^  over 
so  large  a  range  of  magnitudes  as  is  shown  in  figure  6,  although  the  evidence  is  considerable  that  this 
approach  is  valid  in  the  continuum  absorption  wavelengths  ranging  at  least  from  8  to  about  20  pm. 

4.  CONCLUSIONS 


The  work  discussed  has  shown  that  anomalous  electromagnetic  absorption,  at  least  in  the 
IR  continuum  absorption  wavelength  region,  can  be  quantitatively  described  as  being  due  to  distributions 
of  neutral  molecular  clusters  in  water  vapor  and  moist  air  whose  populations,  sizes,  dependencies  on 
humidity  and  temperature,  and  other  parameters,  also  have  been  modeled  here.  The  ratio  modeling 
equations  (Eqs  12-17  and  19-21)  produce  results  in  excellent  agreement  with  observation  and  show 
the  complete  range  of  variation  in  response  to  s  and  0,  not  just  the  absorption  at  constant  partial 
pressure  as  approximated,  for  example,  by  Eq  18.  These  new  equations  are  based  on  the  physics  of 
water  vapor  and,  therefore,  they  should  permit  improved  modeling  of  atmospheric  IR  transmission. 
They  are  recommended  for  evaluation  by  experimentalists  and  for  inclusion  in  existing  atmospheric 
transmission  models  as  warranted.  Since  the  cluster  species  believed  responsible  for  anomalous 
electromagnetic  absorption  are  well  described,  it  should  be  possible  for  workers  in  other  wavelength 
regions  than  the  IR,  e.g.,  in  millimeter  waves  and  microwaves,  to  compute  and  to  experimentally 
verify  absorption  attributable  to  modes  of  neutral  clusters  active  in  the  absorption  at  those 
wavelengths. 
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APPENDIX 

MICROPHYSICS  OF  CLOUDS  AND  NUCLEATION  THEORY 


The  following  discussion  is  not  essential  to  the  development  of  the  equations  given  in  the 
text,  which  are  obtained  directly  from  experimental  data  discussed  there.  This  appendix  is  included 
for  readers  interested  in  the  microphysics  of  clouds  and  nucleation  theory  (gas  to  particle  conversion), 
since  the  interpretations  presented  here  can  be  supported  by  the  various  experimental  data  but  are 
decidedly  contrary  to  existing  theories.  The  agreement  between  these  interpretations  and  observation 
remind  us  that  any  theory  is  only  as  good  as  its  agreement  with  the  data.  This  includes  traditional 
theories,  however  ingrained  in  our  thinking  they  may  have  become  with  the  passage  of  time.  The 
reader  is  encouraged  to  consider  these  arguments  open-mindedly.  If  large  populations  of  neutral 
molecular  clusters  exist  in  water  vapor,  they  surely  must  serve  as  condensation  sites  for  the  growth  of 
droplets  under  suitable  conditions.  Thus  there  would  be  no  need  to  explain  these  sites  as  originating 
from  monomer  collisions  according  to  classical  kinetic  theory,  which  in  any  case  is  imprecise  for  the 
vapors  of  hydrogen-bonded  substances  like  water.  In  this  section,  an  alternative  explanation  is  given 
which  fits  all  observations,  but  for  which  no  further  justification  is  claimed. 

Mass  spectra  and  other  measurements  require  that  the  possibility  be  considered  that,  at  the 
boiling  point,  critical  conditions  are  attained  and  nucleation  and  recondensation  occur  continuously 
between  the  vapor  and  liquid  phases  of  water.  If  essentially  all  vapor  is  clustered  at  the  vapor/liquid 
interface  at  the  boiling  point  (s  =  1.0,  100°C,  1  atm),  then  the  maximum  partial  pressure  fraction  for 
freshly  evaporated  water  vapor  at  the  liquid  interface  at  1  atm  under  other  conditions  than  the  boiling 
point  is: 


(sXpV 

(n_)v  * - 

cv  (1X760) 

where  s  is  the  saturation  ratio  and  p°  is  the  saturation  vapor  pressure  at  temperature  8. 


(A-l) 


This  simple  equation,  with  standard  vapor  pressure  tables,  allows  the  construction  of  what 
might  be  called  equilibrium  phase  diagram  for  vapor-phase  clusters  in  contact  with  liquid  water,  as 
drawn  in  figure  A-l,  for  1  atm  total  pressure.  The  values  of  (nc)v  given  by  Eq  A-l  and  shown  in 
figure  A-l,  are  found  to  be  in  good  agreement  with  estimates  of  the  cluster  fraction  from  measurements 
discussed  in  the  text,  provided  that  liquid  water  is  present  in  the  system.  Since  by  this  interpretation 
all  excess  vapor  must  be  condensed  at  s  =  1.0,  a  vertical  condensation  line  is  shown  there.  Figure  A-l 
shows  that  a  surprisingly  large  fraction  of  water  vapor  is  clustered  even  at  temperatures  and  humidities 
well  below  the  boiling  point,  provided  that  liquid  water  is  present  in  the  system  and  is  in  equilibrium 
with  the  vapor.  For  example,  one  percent  of  all  vapor  can  be  clustered  at  40°C  for  a  saturation  ratio 
s  =  0.15  (15%  RH).  But  this  entire  fraction,  as  discussed  in  the  text,  is  not  effective  in  absorption 
except  near  the  peak  of  the  IR  continuum  absorption  spectrum  (30-50  Jim).  Therefore  measurements 
of  electromagnetic  radiation  absorption  by  neutral  water  clusters  in  the  atmosphere  frequently  infer 
that  cluster  populations,  and  hence  the  fraction  of  vapor  clustered  (nc)v,  are  smaller  than  actually  is 
the  case  and  can  be  demonstrated,  e.g.,  by  measurements  of  the  dissociative  ions  of  the  neutral 
clusters  by  electrical  conductivity.  The  diagonal  line  with  points  in  figure  A-l  and  the  dashed  nearly 
vertical  line  now  will  be  considered. 
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Phase  Diagram  of  Neutral  Clusters  Constructed  from  Eqs  A-l  and  A-2  and  Standard  Tables 
of  Water  Saturation  Vapor  Pressure  According  to  the  Methodology  in  the  Appendix 


If  nucleation  of  neutral  clusters  in  the  vapor  can  occur  at  s  *  1  (an  assumption  long  made 
by  atmospheric  cloud  macrophysicists,  although  without  theoretical  justification),  then  so-called 
supersaturation  need  not  exist  to  explain  nucleation,  as  all  nucleation  phenomena  could  be  explained 
as  temperature  rather  than  partial  pressure  phenomena.  Wilson,1  *2  had  he  known  about  hydrogen 
bonding  and  the  existence  of  equilibrium  neutral  cluster  populations  in  water  vapor  when  his  classical 
adiabatic  expansion  (cooling)  cloud  chamber  experiments  were  performed,  could  have  explained  his 
results  as  due  to  changes  in  temperature  alone.  Saturation  vapor  pressure  is  a  fundamental  physical 
constant  of  water  at  any  temperature,  and  the  concept  of  exceeding  it  arbitrarily  to  explain  new 
results  by  “supersaturations”  of  four  or  five  must  have  caused  concern.  Indeed,  Wilson  had  misgivings 
about  his  observations.3  While  many  empirical  modifications  to  nucleation  models  and  equations  have 
proved  necessary  over  the  years  to  secure  good  agreement  between  theory  and  cloud  chamber  obser¬ 
vation,  these  modifications  have  served  only  to  secure  this  agreement;  they  do  not  prove  that  the 
supersaturation  concept  is  a  correct  one. 

The  interpretation  of  the  boiling  point  as  the  critical  saturation  point  leads  to  simple 
equations  that  are  found  to  fit  experimental  data  more  closely  than  do  traditional  models  based  on 
supersaturation.  For  example,  if  one  considers  the  temperature  drop  that  is  required  to  produce  a 
certain  excess  of  water  vapor  that  is  then  available  to  condense  on  clusters  to  grow  droplets  large 
enough  for  optical  detection,4  then  the  “barrier”  conditions  that  must  be  attained  for  droplet 
formation  can  be  given  as  a  simple  function  for  subsaturated  vapor  referenced  to  the  boiling  point: 

/  373\5^ 

‘Barrier”  Conditions:  ln(s)|,ar  =  1  -  I  )  (A- 2) 


Eq  A-2  is  not  unlike  classical  nucleation  equations,  but  it  gives  exact  solutions  rather  than 
approximate  ones  as  in  the  equations  of  J.  G.  Wilson.5  Eq  A-2  is  the  diagonal  line  with  points 
plotted  in  figure  A-l,  and  its  equivalent  curve  for  Ncc  corresponding  to  (nc)v  is  shown  as  a  nearly 
vertical  function  in  figure  3  of  the  main  body  of  this  report.  Eq  A-2  also  defines  the  right-hand 
diagonal  of  figure  4,  in  equivalent  units  for  (C£).  The  dashed,  nearly  vertical  line  in  figure  A-l  is 
this  barrier  equation  corrected  for  Icc  by  using  the  temperature  dependence  of  the  ion  product 
(Eqs  6,  23).  When  considering  the  neutral  cluster  fraction  in  water  vapor  or  the  cluster  populations, 
Ncc,  the  solid  barrier  curve  with  points  in  figure  A-l,  therefore,  is  applicable. 


Perhaps  fortuitously,  it  is  found  that  (nc)v  can  be  calculated  almost  precisely  at  the  barrier 

curve  by  : 


which  requires  that: 


(A-3) 


(A-4) 
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And  thus  if  one  wishes  to  take  the  ratio  of  the  fraction  clustered  at  two  temperatures  other  than 
373°K,  he  obtains: 


In  figure  A- 2,  Eq  9  is  plotted  to  show  the  mean  size  of  the  neutral  cluster  distribution  vs 
s  and  0  from  the  mass  spectral  data  of  Ref.  3  and  the  interpretations  of  the  text.  The  barrier  curve 
from  Eq  A- 2  is  included. 

All  of  the  examples  of  nonlinear  data  discussed  in  the  text  can  be  explained  by  the  barrier 
curve  of  Eq  A- 2.  IR  absorption  vs  s  becomes  nonlinear  as  the  barrier  curve  is  approached  due  to 
nonlinear  changes  in  the  cluster  population.  Thus  modelers  must  be  careful  to  avoid  saturation  ratios 
approaching  the  barrier  in  their  data  collection,  and  spectroscopists  could  confuse  barrier  conditions 
as  due  to  droplet  growth  on  large,  conventional  condensation  nuclei  when  in  fact  it  was  the  cluster 
population  that  was  responding  to  the  barrier.  This  would  be  the  true  beginning  of  what  might  be 
called  droplet  growth,  but  it  is  distinctly  different  from  that  on  other  nuclei.  It  is  the  beginning  of 
accelerated  growth  of  homogeneous  neutral  clusters  already  present  in  the  vapor,  and  it  will  occur 
regardless  of  whether  other,  larger  condensation  nuclei  are  present  in  the  system  or  not,  provided 
that  liquid  water  is  present.  In  steam,  virtually  all  water  vapor  can  be  clustered,  thereby  explaining 
the  enormous  IR  absorptions  and  emissions  reported  by  many  workers.6'8  after  monomer  and  droplet 
corrections  are  made.  Observations  since  Faraday’s  time  have  shown  that  saturated  steam  contains 
enormous  populations  of  (dissociative)  ions,  but  that  dry  superheated  steam  has  far  fewer  ions  present 
If  temperatures  above  100°C  are  to  be  investigated  for  cluster  activity  in  water  vapor,  the  total 
pressure  must  also  be  raised  to  maintain  equilibrium  conditions.  Since  the  boiling  point  keys  this 
activity,  one  might  expect  that  water  vapor  samples  run  in  partial  vacuum  will  give  variable  results 
depending  on  cluster  species  and  populations  present.  The  dissociative  ion  populations  (Icc)  closely 
track  the  neutral  cluster  populations  (Ncc)  in  experimental  measurements. 

Conditions  to  the  left  of  the  diagonal  barrier  with  points  in  figure  A-l  are  difficult  to 
reproduce  in  the  atmosphere  because,  if  liquid  water  is  present  in  the  system,  s  will  increase  to  the 
barrier  while,  if  liquid  is  not  present,  the  neutral  cluster  population  will  decay  to  lower,  nonequilibrium 
(nc)v  s.  A  common  atmospheric  occurrence  is  that  of  approximately  constant  partial  pressure  with 
temperatures  changing  over  a  period  of  hours,  leading  to  changes  in  the  saturation  ratio  (relative  humidity). 
Constant  partial  pressure  is  shown  by  the  dashed  horizontal  line  in  figure  A-2,  and  results  in  very  small 
changes  in  (nc)v.  Measurements  of  the  ion  population  of  atmospheric  air  show  much  lower  Icc  values  in 
aged  moist  air  than  in  air  recently  in  contact  with  liquid  water. 


Figure  A-2.  Mean  Size  of  the  Neutral  Cluster  Population  from  Eq  9 
Derived  from  Mass  Spectral  Data  of  Ref.  3 

These  are  mean  sizes  only,  and  the  full  distributions  include  clusters  of 
considerably  smaller  and  larger  numbers  of  monomers  (c’s).  Sizes  are 
estimated  by  adding  1  to  the  ion  cluster  measured  sizes,  per  the  discussion 
of  Eq  3.  The  barrier  equation  curve  (Eq  A-2)  is  included  for  reference. 
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OFFICE  OF  THE  SURGEON  GENERAL 
Commander 

US  Army  Medical  Bioengineering 

Research  &  Development  Laboratory 
ATTN:  SGRD-UBD-AL  1 

Fort  Detrlck,  Bldg  568 
Frederick,  MD  21701 

Headquarters 

US  Army  Medical  Research  and 
Development  Command 

ATTN:  SGRD-PL  1 

Fort  Detrlck,  MD  21701 

Commander 

USA  Medical  Research  Institute  of 
Chemical  Defense 

ATTN:  SGRD-UV-L  1 

Aberdeen  Proving  Ground,  MD  21010 

US  ARMY  HEALTH  SERVICE  COMMAND 

Superintendent 

Academy  of  Health  Sciences 

US  Army 

ATTN:  HSA-CDH  1 

ATTN:  HSA-IPM  1 

Fort  Sam  Houston,  TX  78234 


US  ARMY  MATERIEL  DEVELOPMENT  AND 
READINESS  COMMAND 

Commander 

US  Army  Materiel  Development  and 
Readiness  Command 

ATTN:  DRCLDC  1 

ATTN:  DRCSF-P  1 

9001  Elsenhower  Ave 
Alexandria,  VA  22333 

Dt rector 

Human  Engineering  Laboratory 

ATTN:  DRXHE-SP  (CB  Defense  Team)  1 

Aberdeen  Proving  Ground,  MD  21009 

Commander 

US  Army  Foreign  Science  l  Technology 
Center 

ATTN:  DRXST-MT3  1 

220  Seventh  St.,  NE 
Charlottesville,  VA  22901 

Director 

US  Army  Materiel  Systems  Analysis  Activity 
ATTN:  DRXSY-MP  1 

ATTN:  0RXSY-CA  (Mr.  Meti)  2 

Aberdeen  Proving  Ground,  MD  21009 

Commander 

US  Army  Missile  Command 
Redstone  Scientific  Information  Center 
ATTN:  DRSMI-RPR  (Documents)  1 

Redstone  Arsenal,  AL  39809 

D I  rector 

DARC0M  Field  Safety  Activity 

ATTN:  DRX0S-C  1 

Charlestown,  IN  47111 

Commander 

US  Army  Natick  Research  and 
Devalopment  Command 


ATTN:  ORDNA-O  1 
ATTN:  DRDNA-VC  1 
ATTN:  DRDNA-VCC  1 
ATTN:  DRDNA-VM  1 
ATTN:  DRDNA-VR  1 
ATTN:  DRDNA-VT  1 


Natick,  MA  01760 
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US  ARMY  ARMAMENT  RESEARCH  AND 
DEVELOPMENT  COMMAND 

Commander 


US  Army  Armament  Research  and 
Development  Command 


ATTN: 

DRDAR-LCA-L 

1 

ATTN: 

0RDAR-LCU-CE 

1 

ATTN: 

0RDAR-PMA  (G.R.  Sacco) 

1 

ATTN: 

DRDAR-SCA-W 

1 

ATTN: 

DRDAR-TSS 

5 

ATTN: 

DRCPM-CAWS-AM 

1 

ATTN: 

DRCPM-CAWS-S 1 

1 

Dover, 

NJ  07801 

US  ARMY  ARMAMENT  MATERIEL  READINESS 
COMMAND 

Commander 

US  Army  Armament  Materiel  Readiness  Command 
ATTN:  DRSAR-ASN  1 

ATTN:  ORSAR-SF  I 

ATTN:  DRSAR-SR  1 

Rock  Island,  IL  61299 

Commander 
USA  ARRCOM 

ATTN:  SARTE  1 

Aberdeen  Proving  Ground,  MD  21010 

Commander 

US  Army  Dugway  Proving  Ground 

ATTN:  Technical  Library  (Docu  Sect)  1 

Dugway,  UT  84022 

US  ARMY  TRAINING  4  DOCTRINE  COMMAND 
Commandant 

US  Army  Infantry  School 

ATTN:  NBC  Division  1 

Fort  Banning,  GA  31905 

Commandant 

US  Army  Missile  &  Munitions  Center 
and  School 

ATTN:  ATSK-DT-MU-EOD  1 

Restone  Arsenal,  AL  35809 

Commandant 

USAMP4CS/TC4FM 

ATTN:  ATZN-CM-CDM  1 

Fort  McClellan,  AL  36205 


US  Army  Chemical  School 

ATTN:  ATZN-CM-AL  1 

Fort  McClellan,  AL  36205 

Commander 

US  Army  Infantry  Center 

ATTN:  ATSH-CD-MS-C  1 

Fort  Bennlng,  GA  31905 

Commander 

US  Army  Infantry  Center 
Directorate  of  Plans  4  Training 
ATTN:  ATZB-DPT-P0-NBC  1 

Fort  Banning,  GA  31905 

Commander 

USA  Training  and  Doctrine  Command 

ATTN:  ATCD-Z  1 

Fort  Monroe,  VA  23651 

Commander 

USA  Combined  Arms  Center  and 
Fort  Leavenworth 

ATTN:  ATZL-CA-COG  1 

ATTN:  ATZL-CAM-IM  1 

Fort  Leavenworth,  KS  66027 

Commander 

US  Army  TRADOC  System  Analysis  Activity 
ATTN:  ATAA-SL  1 

White  Sands  Missile  Range,  NM  88002 

US  ARMY  TEST  4  EVALUATION  COMMAND 

Commander 

US  Army  Test  4  Evaluation  Command 

ATTN:  DRSTE-CT-T  1 

Aberdeen  Proving  Ground,  M0  21005 

DEPARTMENT  OF  THE  NAVY 

Chief  of  Naval  Research 

ATTN:  Code  443  1 

800  N.  Quincy  Street 
Arlington,  VA  22217 

Commander 

Naval  Explosive  Ordnance  Disposal  Facility 
ATTN:  Army  Chemical  Officer  (Code  AC-3)  1 

Indian  Head,  MD  20640 
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Command er 

Naval  Surface  Weapons  Center 

Code  G51  t 

Dahlgren,  VA  22448 

Chief,  Bureau  of  Medicine  &  Surgery 
Department  of  the  Navy 

ATTN:  MED  3C33  1 

Washington,  DC  20372 

Commander 

Naval  Weapons  Center 

ATTN:  Technical  Library  (Code  343)  1 

China  Lake,  CA  93555 

US  MARINE  CORPS 

Director,  Development  Center 
Marine  Corps  Development  and 
Education  Command 

ATTN:  Fire  Power  Division  t 

Quantlco,  VA  22134 

DEPARTMENT  OF  THE  AIR  FORCE 

HQ  Foreign  Technology  Division  (AFSC) 

ATTN:  TQTR  1 

Wrlght-Patterson  AFB,  OH  45433 


Commander 

Aeronautical  Systems  Division 

ATTN:  ASD/AESD  1 

Wr I gh t-Patterson  AFB,  OH  45433 

HQ  AFLC/LOWMM  1 

Wr Igh t-Patterson  AFB,  OH  45433 

HQ.  AFSC/SDNE  1 

Andrews  AFB,  MD  20334 

HQ  AMD/RD 

ATTN:  Chemical  Defense  OPR  1 

Brooks  AFB,  TX  78235 

NOR AO  Combat  Operations  Center 

ATTN:  DOUN  1 

Cheyenne  Mtn  Complex,  CO  80914 

AFAMRL/HE 

ATTN:  Dr.  Clyde  Rep  toggle  1 

Wrlght-Patterson  AFB,  OH  45433 

HO  AFTEC/SGB  1 

Klrtlend  AFB,  NM  87117 


USAF  TAWC/THLO  1 

Eglln  AFB,  FL  32542 

USAFSAM/HR 

ATTN:  Dr.  Robert  Reyes  1 

Brooks  AFB,  TX  78235 

USAFSAM/VNT 

ATTN:  Dr.  F.  Wesley  Baumgardner  1 

Brooks  AFB,  TX  78235 

AFAMRL/HE 

ATTN:  COL  Dan  Johnson  1 

Wrlght-Patterson  AFB,  OH  45433 

AMD/ROE 

ATTN:  LTC  T.  Klngery  1 

Brooks  AFB,  TX  78235 


OUTSIDE  AGENCIES 

Battalia,  Columbus  Laboratories 
ATTN:  TACTEC  1 

505  King  Avenue 
Columbus,  OH  43201 

Toxicology  Information  Center,  WG  1008 
National  Research  Council  1 

2101  Constitution  Ave. ,  NW 
Washington,  DC  20418 

US  Public  Health  Service 

Canter  for  Disease  Control 

ATTN:  Lewis  Webb,  Jr.  1 

Building  4,  Room  232 

Atlanta,  GA  30333 
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